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The eigenvalues of quantum chaotic systems have been conjectured to foUow, in the large energy 
limit, the statistical distribution of eigenvalues of random ensembles of matrices of size A'^ — >■ oo. 
Here we provide semiclassical arguments that extend the validity of this correspondence to finite 
energies. We conjecture that the spectrum of a generic fully chaotic system without time-reversal 
symmetry has, around some large but finite energy E, the same statistical properties as the Circular 
, Unitary Ensemble of random matrices of dimension NaS = Tu/y/2Ad\, where th is Heisenberg time 

fS| ' and y/di is a characteristic classical time, both evaluated at energy E. A corresponding conjecture 

^ ' is also made for chaotic maps. 
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I. INTRODUCTION 



The study of connections between Random Matrix Theory (RMT) and properties of chaotic systems has now a 
' long history [H, A central point is the Bohigas-Giannoni-Schmit conjecture [3]. For systems without time reversal 
P5 ■ invariance, the conjecture states that in the asymptotic high-energy limit the fluctuation properties of eigenvalues 
of quantum chaotic systems coincide with those of the Gaussian unitary ensemble (GUE) of random matrices of 
dimensionality iV — > oo. In particular, the normalized pair correlation function of the asymptotic eigenvalues of a 
chaotic system is conjectured to be 

>' i?2(.)=i-f^V> (1) 



o 

where x is the unfolded distance between eigenvalues (i.e. the mean spacing has been set to one). 
' [ Nowadays, RMT has a large number of applications ranging from disordered systems in quantum mechanics, 

04 to growth processes in statistical physics and zeros of Riemann Zeta function in number theory. It is often used to 
obtain deep conjectures or to compute particular properties of physical systems which are usually beyond the reach 
of exact methods. The validity of these predictions is then tested via experiments or via numerical calculations. 

One of the difficulties concerning such an approach is that often real systems are not in the asymptotic regime 
required by the conjecture, and it may be experimentally quite unnatural to force them into it. Another point is 
related to the numerical tests of the predictions, which due to their asymptotic nature often require difficult large 
scale computations. 

' A relevant question concerns therefore the statistical properties of the finite-energy eigenvalues of quantum chaotic 
systems, and its eventual connections to RMT. In the following we will argue that such a connection exists. More 
precisely, based on a semiclassical analysis of the two-point correlation function we conjecture that the spectrum of 
a generic fully chaotic dynamical system without time-reversal symmetry has, around some finite energy E, the same 
statistical properties of the eigenvalues of the Circular Unitary Ensemble (CUE) of random matrices of dimension 



N,s = , (2) 



where 



di = — lim 

2 ri— >oc 



Tp<Ti 



r2|det(Mp-/)| 



(3) 



In Eq. (HI), Th = h'p{E) is the Heisenberg time, where p(i?) is the smooth part of the quantum density of states. The 
factor \/d^ is a system-dependent classical time expressed in terms of a sum over all the classical periodic orbits p of 
the system at energy E, Tp is the period of the periodic orbit, and Mp it's monodromy matrix (/ is the identity matrix). 
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The sum is over all periodic orbits of period less than ri , including the repetitions rp when orbit p is non-primitive 
[vp = 1 for primitive periodic orbits). Although not explicitly indicated, N^f should be interpreted as the closest 
integer approximation of the r.h.s.. In the semiclassical limit E oo, A'eff oo, and we recover the asymptotic CUE 
result (that coincides with the Gaussian Unitary Ensemble), e.g. we recover the Bohigas-Giannoni-Schmit conjecture. 

A similar problem may be considered for another type of systems, when the classical chaotic dynamics is described 
by a canonical map. These systems have been quite useful as schematic models for studying the quantum-classical 
correspondence, since a chaotic dynamics is possible in one dimension. When their phase-space is compact, they also 
allow to avoid truncations of the quantum mechanical model, since their evolution operator is given by a unitary 
matrix of finite dimension N, where the limit N ^ oo describes the semiclassical limit. Connections between chaotic 
maps and random matrix ensembles have been already established [5| in this limit. Beyond their formal character, 
maps have also been quite useful to analyze real systems since, sometimes, they provide quite good approximations 
to their dynamics @ . Well known examples of such maps are the standard map T| , the baker map 8] , or the kicked 
Harper map 0. 

We have also made a semiclassical analysis of the statistical properties of the eigenphases of quantum chaotic maps 
without time-reversal invariance. For models described by a finite unitary matrix U of size N, we conjecture that the 
spectrum of these chaotic maps have the same statistical properties of CUE random matrices of dimension 

N 

NeS - , , (4) 

VI + 24ci ^ ' 

where 

1 

Z n* — ^oo ^ — ^ 
n—1 

is a finite real number that depends on the system considered (a„ = trt/" are the traces of powers of the evolution 
operator U). 

The present work may be viewed as a continuation of the results obtained in Ref. [lo| , where it was conjectured that 
zeros of the Riemann zeta function located around z = 1/2 + IE' have the same statistical properties as the eigenvalues 
of CUE random matrices of effective dimension 

log(i;/2^) 

TVcff = (6) 



where /3 = 1.57314, is a constant depending on sums over prime numbers (see the discussion in Section |V|) . 

The paper is organized as follows. In Section|lI]general expressions are given for the density of states, its semiclassical 
approximation, and the two-point correlation function. For chaotic Hamiltonian systems, the method developed by 
Bogomolny and Keating [l^ is used to express the two-point correlation function in terms of a classical function, 
which is formally exact. A similar method for maps is also developed. 

In Section IlIII the classical functions obtained in Section |ll] are expanded for small arguments, which corresponds 
to the semiclassical limit, keeping all orders. 

In Section HVl we first consider the exact two-point correlation function of the Circular Unitary Ensemble at finite 
dimension TV. Secondly, the expansion of the classical function from previous section is used to compute the two-point 
correlation function up to second power in Heisenberg time for chaotic Hamiltonian systems without time reversal 
symmetry. Analogously we obtain a result for maps up to order N'^. Comparing these results to those of CUE, we 
derive the effective dimension. 

Finally, in Section|V]we check the results for the Riemann Zeta function, making some general concluding remarks. 



II. SEMICLASSICAL CALCULATION OF THE TWO-POINT CORRELATION FUNCTION 



In the semiclassical limit, the density of states 

p{E) =Y,S{E - E,) (7) 

may be split into a smooth plus an oscillatory part p = ]) + p. The smooth part admits an expansion in powers of 
energy called the Weyl series, that depends on global properties of the system such as volume, surface and topology 
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of the classical phase space, disregarding the integrable or chaotic nature of the dynamics. The oscillatory part may, 
in contrast, be written 

p{E) = 2 ME) cos [Sp{E)/h + lyp] , (8) 
p 

where as in Eq. ([3]) the sum is over all classical periodic orbits p of the system at energy E, Sp is the corresponding 
action, and i^p is a phase called the Maslov index that depends on the number of turns of the stable and unstable 
varieties along the periodic orbits. 

The value of the amplitude Ap depends on the stability of the periodic orbit. If the orbit is unstable, the amplitude 
is given by the Gutzwillcr trace formula 11] 

A = _ (9) 

^ 27r?irp^|det(Mp-/)| ' 

Several statistical spectral functions are of interest. One of the most basic ones, on which we will concentrate on, 
is the two-point correlation function 

7^2(ei, £2) = {p{E + €i)p{E + £2)) . (10) 

The brackets indicate an energy average defined as 



ifiE)) - -rr; / f{E')dE' , (11) 



where the energy window AE is small compared to E but it is however large enough to be statistically meaningful 
(usually AE > 1/p). 

A. Hamiltonian systems 

Semiclassically, the r2 function may be computed by introducing the previous semiclassical approximation into the 
definition (|10p . The energy average of the double sum over periodic orbits splits into a diagonal and an off-diagonal 
part. By the Hannay-Ozorio de Almeida sum rule the diagonal part can be easily estimated. The main problem is to 
compute the off-diagonal term. Actions of different periodic orbits have correlations that are not yet fully understood. 
To overcome this problem we follow the pioneering work of Bogomolny and Keating , where it is argue that the 
two-point correlation function of a chaotic system tends to the random matrix prediction in the semiclassical limit. 
We sketch here some details of the calculation in order to compute next to leading order corrections, taking into 
account all the approximations. 

We begin by defining for chaotic systems a truncated density of states and counting function 

^^h^^^,rp\detiMp~I)\'^' 

_ 2 .5-^ gi(Sp/;i+!^p) 

K*iE) = V 777+c.c. , (13) 

2^\f^^.rp\detiMp^I)\'^' 

If T* is of the order of Heisenberg time, and 

J\fr.{E) ^U{E)+Afr'{E) = n+ 1/2 (14) 

then the energies E = En{T*) obtained by inverting the latter equation with n = 0, 1,2, . . . provide a semiclassical 
approximation of the exact spectrum. Once we have these approximate eigenvalues we define a bootstrapped density 
of states 

Dr^{E) = Y,S[E~En{T*)] (15) 

71 

As r* — )■ oo, it is obvious that Dr*{E) p{E). Rewriting the sum and using Poisson summation formula 
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Dr* {E) = pr- (E) (-1)'' e^P[i k Nr* {E)] . (16) 

fc— — oo 

In contrast to the truncated density of states p^. , this function has the correct delta-peak structure with the same 
periodic orbit information. Introducing this approximation into the definition of the two-point correlation function 
we obtain the following expressions 

r,{e,,e,)=t{E)+rf-^\e,,e,)+r'^f\e,,e2) , (17) 

with 



rf^^\ei,e2) = (p.. (i? + ei)p.. + £2)) , (18) 
ri°''\ei,e2) ^ ^ ^4 fc ' i^M i^nhUr- {E + a) - lii^k^Mr' {E + 62)]) ■ (19) 

Using the definition (fT2|). we obtain a double sum over periodic orbits. Expanding the action Sp{E+ej) = Sp{E)+Tpej, 
and considering that periods up to t* are uncorrelated, it follows that 

(diag). ^ 1 'f- ffpT^ COs(Tp(ei -62)/^) 

^(^^-^^)-^E^ rg|det(M,-/)| ■ 

We remark that the sum is made over all periodic orbits (primitive, repetitions, symmetric partners). The factor gp 
is usually one, except when the orbit p has a time reversal symmetric partner having the same period and action, in 
which case gp = 2. 

For the off-diagonal term (|19p . after expanding the counting function in the exponential, the average over energy 
kills all terms with fci ^ k2, simplifying to 

with 



(ei,e2) = (exp i27rA:[A/;. (£: + ei) - A/;, (i; + £2)] ) 



(22) 



It is known that for chaotic systems, the fiuctuating part of the counting function is gaussian distributed [13[. Thus, 
we can use the relation (exp[i/(i?)]) = exp[— (/^(£^))/2], valid for any Gaussian random function f{E) with zero 
mean. Finally, applying again the diagonal approximation for the correlator of the counting function 



$fc(ei,e2) 



A^. ((61 - 62)/^) 

A..(0) 



(23) 



where 



A,-, (u) = exp 



Tp<T* 



gp cos (rpu) 
r2 |det(Mp 



(24) 



is a classical function containing all dynamical information of the system. The diagonal term can be also written 
in terms of this classical function, and noting that both diagonal and off-diagonal terms depend only through the 



difference e = ei — £2, we finally rewrite the two-point correlation function as r2(e) = 



^(diag) 



with 
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^2 (e) = - 



E 

k=l 



cos(27rA;ep) 

fc2 



A,, (e/fi) 



A,.(0) 



(25) 
(26) 



In order to compare the correlation functions of different systems, or at different energies, to random matrix tlieory 
results, it is necessary to unfold the spectrum. The natural local energy scale is the mean level spacing S = 1/p, and 
we define the unfolded two-point correlation function as 



^2(2;) = ir2 (3 

where a; is a dimensionless number. Similarly, R2 can be split into a diagonal and off-diagonal part 



(27) 



R2ix) = l + Rf^^\x) + Rr'{x) , 



(off). 



(28) 



with 



1 



Rr\x) 



47r2 
' 27r2 ^ 

k=l 



logA^. (x/hp) 



cos{2nkx) 

fc2 



A, {.r/m 



A.*(0) 



(29) 

(30) 



B. Maps 

The quantization of classical canonical maps gives a unitary evolution operator U whose eigenvalues e'^^ are on 
the unit circle. For maps in a compact phase space, the evolution operator is represented by a finite N dimensional 
unitary matrix, where N = S/2iTh, and S is the phase-space area. We define the density of states over the unit circle 
as 



00 N 

'^W= E Y.^{0-e,-2T,k) , (31) 

where the eigenphases 6j are for convenience restricted to the basic period, < < 27r. By definition p(9) is a 2tt 
periodic function, p{9 + 27r) — p{9). As for continuous systems, the density of states can be split into a smooth and 
fluctuating part p{9) = p{9) + p{9). The smooth part is uniform over the unit circle, 'p{9) = N/2tt. The fluctuating 
part can be written as 



00 

^(^) = ^E[<e'"'+a„e-™^] , (32) 



27r ^ 

n=l 



where a„ = trC/" arc the traces of powers of the evolution operator. By unitarity of U we have a_„ = tiU " = a* . 
Analogously, a cumulative density of states, or counting function, can be defined as 

A/-(0) = / p{9')d9' =77(9) +M{9) . (33) 
^0 



The smooth part is 



J7{0) = ^(^ + C , (34) 
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where 



and the fluctuating part is 



C 



-E 



Im(a„) 



^ OO 



Both the counting function and its smooth part satisfy the translation relation M{d + in] 
fluctuating part is strictly 2it periodic. 

As for hamiltonian systems, the two-point correlation function is defined as 



(35) 
(36) 

= U{e) + N, while the 



^2(ei,e2) = {p{e + e^)p{e + e2)) , 
where the brackets now indicate an average over the phase 



(37) 



In terms of the traces of ?7", the average for r2 can be computed analytically 



(38) 



'^2(ei,e2) 



^ OO 

^E 



47r2 



in(ci-C2)l 



(39) 



n=l 



It has also been conjectured that for fully chaotic classical maps, in the semiclassical limit, the eigenvalue statistics are 
the same as for circular ensembles of random matrices. In order to see this semiclassically, it is necessary to estimate 
the behavior of the traces a„. Tabor [31 developed for chaotic maps a trace formula similar to Gutzwiller's, 



trC/" = 



E 



ne 



i27rrp(JVSp-i/p/4) 



rp|det(A/p-/)|i/2 



(40) 



where the sum is made over all periodic orbits p of period n. To compute |a„p, necessary for the two-point correlation 
function, a double sum on periodic orbits is needed. For periods rt 3> 1 the exponential proliferation of periodic orbits 
in chaotic maps makes that terms in the double sum coming from different orbits approximately cancel, provided 
that actions are uncorrelated, and only the diagonal term survives. In the regime 1 n N, the values of these 
traces are independent of N. For small n, the few off-diagonal terms in the double sum on periodic orbits still depend 
on N, but averaging over an ensemble of different iV values around a large iVo they will also cancel. It follows that, 
semiclassically. 



where gp — 2 when the orbit p has a symmetric partner, otherwise gp = 1. By a classical sum rule for chaotic maps, 
the sum can be further simplified, giving 



|a„P«gn. (42) 

However, it is known that for periods n of order N, non-trivial correlations among actions of different periodic orbits 
become important and the diagonal approximation starts to fail. To overcome this problem we adapt the method 
developed by Bogomolny-Keating to the case of maps, first defining a truncated density of states and counting function 
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N 

PM{e) = 7r + pM{s) , 

N 

^M{e) = Tr^ + c + NM{e) , 

zvr 



(43) 
(44) 



with 



M 

ZTT — ' 



M 



27ri ■^-^ 



'•n An9 



(45) 
(46) 



The truncation number M should be of order N, in which case the exact eigenphases 0j can be well approximated by 
the condition 



J\fM{0j{M))= J + 1/2 



(47) 



for I < j < N and < Oj{M) < 2tt. With these approximated eigenphases we define a bootstrapped density of states 



oo N 

Dm{0) = Y^i.^- ^jiM) - 27rfc) 

A;— — OO j—1 

Using definition ()47p and Poisson summation formula, it can be expressed as 



(48) 



Dm{0) - pm{9) V (-l)'=e'2-^-('') 



k— — oo 



(49) 



Similarly to Hamiltonian systems, this function has the correct delta-peak structure in contrast to the truncated 
density of states . Introducing this approximation into the definition of the two-point correlation function we have 

r2(ei,e2) = {DM{e + ei)DM{e + e2)} 

= {pm{0 + ei)pM{e + 62) J2 (-1)'^'"''' exp[ i27:hAfMiO + d) - i27rfc2A/'M(0 + £2)]) 



ki ,k2 



. (diag) / N (off)/ N 

■^+^2 (ei,e2)+r-^ (£1,62) 



(50) 



In the diagonal part the average on can be done analytically, giving an equation similar to (j39p . although the sum 
is truncated at M and can be rewritten as 



(diag), >^ 



1 32 



Att^ deide2 ^ 

n—l 



^(ei-e2) _|_ g-i"(ei-e2)j 



(51) 



For the off-diagonal term 



(off), s 1 92 



^2 (ei,e2) 



J2 M:4^(exp[ i27:hMM{0 + ei) - i2nk2MMi0 + £2)]) , 



(52) 



it can be shown that terms with ki 7^ k2 cancel after averaging over 9, giving 
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(off) ( \ 



47r2 Sei^ea A:^ 



(53) 



with 



$fc(ei,e2) = (exp[i27rA;[^M(^ + ei)-Mw(^ + e2)]J; • (54) 
Considering that for chaotic maps the fluctuating part of the counting function is also gaussian distributed, we have 



$fc(ei, ea) « exp -27r^fc' ( {Mm{Q + ei) - AAm(^ + 



For maps the correlator of the counting function can be computed exactly in terms of the traces o„ 



M 



exp 




(A/'M(^ + ei)7VM(^ + e2)) 


exp 


47r2A;2([Mw(^)]2)' 





Tl=l 



Introducing the function 



^ la 12 



2 > — -- cos ne 



n=l 



AM(e) = exp 

both the diagonal and off-diagonal parts of the two-point correlation function are expressed as 

1 52 



^(diag)/ X _ 



„(°ff) 



(e) = - 



(e) 



47r2 ae2 

2 °o 



log Am (e) , 



fc=i 



(0) 



(55) 



(56) 



(57) 



(58) 

(59) 



The unfolding for maps is simpler as compared to generic dynamical systems. The smooth density of states is 
uniform over the unit circle, and one has to divide by a constant. The unfolded variable is x = e/p giving e = 2'kx/N, 
and the above expressions become 



R. 



(diag) 



{x) 



1 52 



logAM(27ra;/7V) 



47r2 9a;2 

2^ T2 cos(27rfca;) 



27r2 dx^ fc2 
fe=i 



Am(0) 



(60) 
(61) 



where the total unfolded two-point correlation function is R2{x) = 1 + r''2^'^^\x) + R^^\x). The latter equations 
show that in the semiclassical limit N ^ 00 the correlation function depends on the behavior of the Am function at 
small arguments. We stress that these expressions are completely analogue to those for Hamiltonian systems. 

In principle. Am is a quantum function due to its dependence on the traces |a„|2. However, as discussed before, 
this function may be semiclassically approximated by 



Am(s) 



A^^)(.)=exp 



M 



EE 



gp cos(ns) 
r2|det(Mp-7)| 



(62) 



which is a purely classical function. 



9 



III. EXPANSION OF THE CLASSICAL A FUNCTION 

From the previous section, for Hamiltonian systems we observe that the two-point correlation function is fuhy 
determined by the classical function A^.. (u). We are interested on the behavior of R2{x) for x of order one. Then the 
argument of the classical A^.. function is 

X 2ttx 

u = — ^ , 63 

hp th 

where Th is the Heisenberg time. In the semiclassical limit this time usually goes to infinity, implying m — ^ 0. For this 
reason we are specially interested in the behavior of A^» for small arguments, and we should expand the function 
around u = 0. In order to obtain this expansion we write the function in an integral form 

logA^.(u) = 2/ \KD{T)cos{uT)dT , (64) 







where 



r2 



is the diagonal form factor. For chaotic systems this sum may be estimated by Hannay-Ozorio de Almeida sum rule, 

KD{T)^gT. (66) 

This approximation is valid for r 3> Tj^in, where Tmin is the period of the shortest periodic orbit. The g factor is 
2 for systems with time- reversal symmetry, otherwise g — 1. For sufficiently large r*, we can consider an arbitrary 
intermediate time ri , such that Tmin <C ti <C r* , then 

/ N „ ■^-^ Qv cos(uT„) „ cos(mt) , , , 

The last integral can be solved in terms of the cosine integral function 

where 7 — 0.5772 ... is the Euler-Mascheroni constant, giving 

Expanding the cosine into the periodic orbit sum for small u, and using the Taylor series for the cosine integral 
function, we obtain 



logA^.(M) = 2g 

with 



do + ^ d/u^' - log u - 7 + Ci(uT*) 



1=1 



(70) 
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lim 



\p 9p/9 
^Z.^^2|det(Mp-/)| 



logTl 



i9p/9)Tp 



E^^^2|det(A/p~/)| 2/ 



(71) 



(72) 



The latter constants depend on classical information. To compute At*(0), necessary for the off-diagonal part of i?2, 
we easily solve the integral of Eq. (p7)) arriving to 



logA,.(0) = 2.g[do+logT* 



(73) 



For maps, in the semiclassical limit N ^ oo, provided that there exist a n* such 1 <^ n* ^ M, the sum over traces 
in Eq. ([57]) can be split 



M 



logAM(e) = 2^cos(ne)^ + 25 ^ 



cos(ne) 



(74) 



n— 1 n— n*+l 

where in the second sum we have used the classical summation rule (j42p . |a„p ~ gn. After rearranging terms we have 



log AM(e) = 6„cos(ne) + 2g 



Ecos(ne) 



(75) 



with 



_ Kr _ 1 



5"' " ^^p|det(Mp-/)| 



(gp/g) _ 1 
n 



(76) 



These coefficients measure the deviation of traces with respect to the classical sum rule. Therefore it is expected that 
6„ — > as n — > cx). Moreover, if 6„ decreases sufficiently fast as n increases, the first sum becomes independent of n*, 
and we can expand the cosine for small e. The first sum in (j75p is 



bn cos(ne) = ^ Qe^' 



(77) 



n=l 



with 



(-1) 
(20! 



lim bnT 



(78) 



n=l 



If M —i' oo, the second sum in ((75|) has an analytic solution. The error due to a finite M can be estimated through 
the Euler-Maclaurin formula 



M , . 

cos(ne) 



= -log|2sin(e/2)|+FM(e) , 



(79) 



with 



FM{e)^- V £^^Ci(M6) + 0(M-i) 



n=M+l 



(80) 
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For fixed e, as M — )■ oo, Fm{£) is small. Finally 



1 



|2sin(e/2)|29 
For e = 0, we evaluate Eq. ([75]) directly 



exp { 2g 



CQ + ^cie"^^ + FM{e) 



1=1 



log Aa,/(0) = 2gJ2 bn + 2gJ2-= 2g[co + Hm] 



n—1 n—1 

where Hm is the harmonic number. It is known that for large M the harmonic number admits the expansion 

Hm = log Af + 7 + -i- - -1- + 0{M-^) , 



and finally 



2M 12M2 



(81) 



(82) 



(83) 



AM(0) = M29exp 25 



1 1 



Co + 7- 



2M 12A/2 



0{M- 



(84) 



IV. COMPARISON TO RANDOM MATRIX ENSEMBLES 

In RMT, correlation functions and other related statistics are computed averaging over the ensemble, instead of 
averaging over the spectrum for a particular realization of the ensemble (although it is known that for sufficiently 
large matrices both results coincide). In several cases, exact analytical expressions are obtained. We are interested 
specially in random matrix ensembles of finite dimension. For the CUE of dimension N, a simple exact expression 
for the two-point correlation function is [l3| 



R2{x) = 1 



sin(7ra;) 



Nsin{TTx/N) 



As iV — > cx) it expands as 



R2ix) = 1 



sin(7ra;) 



sin^(7ra;) 



157V4 



(85) 



(86) 



in terms of even powers of the dimension N. The leading order is the usual RMT result for the R2 function, that also 
coincides with that of the Gaussian Unitary Ensemble (GUE) in the limit N 00. 

A. Hamiltonian systems 

Inserting Eqs. ([70)) and ([75]) for A^. into Eqs. ([2^ and ([50]) we can compute the unfolded two-point correlation 
function in terms of classical quantities di . In order to simplify the equations and compare them to random matrix 
results we have to make further approximations: 

i) For small but fixed u, if we consider r* — >■ 00, the term Ci(ur*) in Eq. (|70l) vanishes. 

ii) As we considered t* ~ Th, and At.(O) ^ (t*)^^, we keep only the term fc = 1 in the sum of the off-diagonal part. 

iii) For the second derivative in the off-diagonal part we neglect the derivatives of A{u)t-' because they give higher 
order terms in Th. 

The two-point correlation function for dynamical systems, up to second order in Heisenberg time, simplifies to 



R. 



(dii 



1 „ , /27r 
^ + 2di — 

X'' \Ta 
.2g 



(87) 



4°''\x) = 2 (2^^) [1 - 2sin2(^x)] exp [2gd,{27rx/T^r + ©(r^^)] , 
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Moreover, following the same argument as in Ref.|T2|, for systems without time-reversal symmetry (g = 1) we can 
take T* = Tn/e'*', in order to exactly cancel the divergence in the diagonal part by the ofF-diagonal part as x — >■ 0, 
arriving to 



R2ix) = 1 



sin(7ra;) 



TTX 



^sin2(^x)+0(r„4) 



(89) 



We see explicitly that in the limit Th — >■ oo we recover the asymptotic CUE or GUE result, and we have a first order 
correction for finite Heisenberg time depending on periodic orbits through the coefficient di. Moreover, comparing to 
Eq. (|86| for the CUE ensemble, we see that the functional form in x of the first correction is the same. Equating the 
constants it follows that a generic chaotic dynamical system without time-reversal symmetry has the same two point 
function as a finite CUE ensemble of dimension 



iVcff = 



(90) 



where the value of di was given in Eq. ([3]). 

An estimate of the characteristic time associated to di may be obtained as follows. The sum over periodic orbits 
in Eq.Q may be approximated by 



drxir) 



|det(M^ - J)| 



(91) 



where x(''") — exp(/it„pT)/r — exp(ht„„T /2) /2t — ... is the smooth density of primitive periodic orbits of period r, and 
/itop is the topological entropy (cf Ref.[15]). Using the approximate sum rule |det(Af7. — /)| w exp(/it„pT), it follows 
from Eq. (pij) and the definition of di that 



di 



1 

top 



(92) 



B. Maps 



We compute the unfolded two-point correlation function for maps inserting Eqs. ([5T|) and ([M)) for Am into Eqs. 
Ml and (pT|). 



9 
"27r2 



R 



(off) 



(x) = 



sin^inx/N) 
2 



2ci 



[2MeTsin(7rx/iV)]29 



[1 — 2sin^(7rx)] exp 



(93) 
(94) 



All the information about the map is contained in the constant ci, that can be computed semiclassically through 
the traces a„. For maps without time- reversal symmetry ((7 = 1), if we take M = N/e^ , and if we first consider the 
constant ci = 0, we arrive to 



i?2(x) 



sin(7ra;) 



7Vsin(7ra::/iV) 



(95) 



which is exactly the same exact result ([85]) for A^-dimensional CUE matrices. This would show that evolution operator 
matrices of chaotic maps without time-reversal symmetry behave as CUE random matrices of the same dimension. 
Nevertheless, considering that constant ci is finite, up to second order in N, the two-point correlation function is 



R2(x) 



sin(7ra;) 



TTX 



^^.in^inx)+OiN^^) 



(96) 
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Comparing to Eq. (l86l) . we conclude that corrections to the asymptotic result are the same as a CUE matrix of 
effective dimension 

The value of constant ci for maps was explicitly given in Eq. ([5|). Similarly to Hamiltonian systems, an estimation 
for Cl can be done from the density of periodic orbits. Here the integral over time is replaced by a sum over integer 
periods n, giving 

C-'Uop/2 

~ 4(1 _ e-''.=p/2)2 • 
where now the topological entropy /itop is a dimensionless number. 



V. CONCLUSIONS 



The results obtained here for the two-point correlation function may, hopefully, be extended to arbitrary n-point 
functions (see the analoguous discussion for the Riemann Zeta function in Refs. [13] and [l^). This extension would 
imply that all statistical properties of eigenvalues of chaotic systems without time reversal symmetry will also agree 
with that of circular random matrices of effective dimension Ncs — T^j \flAd\. In particular, this should be the case 
for the nearest neighbor spacing distribution. 

The example of the Riemann Zeta function allows to make explicit comparisons with our results, due to the huge 
numerical data on zeroes and explicit periodic orbit information [l^. Following [l||, the prime numbers p label the 
primitive periodic orbits of period Tp = logp, with |det(Afp — /)| = ]f . With these replacements, we get 



d\ = — lim 

2 Tl— >OC 



2 ^ 

r logp<Ti 



log P 



(99) 



where the sum is made over primes p and integers r. This limit gives di = 0.78657, which is exactly half the value of 
(5 = 1.57314 . . . obtained in [10|, and gives in Eq.® the correct effective dimension. In Ref.[l3], the nearest neighbor 
spacing distribution was numerically checked for different heights E in the spectrum with high precision. 

For general chaotic dynamical systems, the computation of c?i needs explicit classical information on periodic orbits. 
However, for ti ^ Tmin, the sum will converge towards /2. Therefore, the numerical computation of di relies mainly 
on the computation of the contribution of short periodic orbits. 

For maps, the constant ci can be computed using the quantum values of the traces or, alternatively, their semiclas- 
sical approximations in terms of periodic points. 

One may wonder if CUE asymptotics may be used instead of CUE. Our analysis indicate that the correct ensemble 
is the latter since, due to a non uniform density of eigenvalues, the CUE has a different asymptotics compare to 
Ea.(|551). In particular, the first correction goes like {—1)^/N, which is not derived from our semiclassical expansion. 

It would be also desirable to extend this study to systems with time-reversal symmetry, and compare to the Circular 
Orthogonal Ensemble for finite N, although expressions are more complex even in the limiting case N ^ oo. 
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